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BOUNDARY LAYER ANALYSIS OF A CENTAUR 


STANDARD SHROUD 
by W. R. Hingst and C. E. Towne 
Lewis Research Center 

SUMMARY 

An analytical boundary layer investi^tion was carried out in 
conjunction with an experimental wind tunnel test to determine the 
discharge characteristics of the Centaur shroud ascent vent system 
on the Titan/Centaur launch vehicle. This involved estimating the 
effect of the local boundary layers on the vent discharge. The growth 
of the boundary layer along the vehicle was influenced by the inter- 
action with flanges protruding into the flow over the body and by the 
loi^tudinal corrugations in the vehicle surface. In addition, boundary 
layer calculations were made for a 3 percent modal of the launch vehi- 
cle, and were compared with e:q)eridiental results. The effect of the 
flange interacting with the boundary layer was determined using a con- 
trol volume technique. The calculation of boundary layer growth over 
the corrugations was estimated by using a surface rou^ness in the 
boundary layer calculations that would simulate the effect of the corru- 
gations. Calculations were made for both the vehicle and model for 
Mach numbers ranging from 0. 8 to 1. 56. A comparison of analytical 
results with experimental data shows very good agreement for the 
supersonic Mach numbers. At the lower Mach numbers the agreement 
between experiment and analysis was not as good. 
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INTRODUCTION 

During preparation for the Titan/Centaur-Vikii^ proof fli^t, it 
was necessary to verify the adequacy of the Centaur shroud ascent vent 
system on the Titan/Centaur launch vehicle. The ascent vent system 
was designed to vent the internal air of the shroud during the ascent 
phase of the mission to prevent significant pressure loading on the vehi- 
cle. The problem essentially resolved itself into the ability to estimate 
the correct discharge coefficient for the vents while the vehicle was in 
the transonic flight regime. This involved estimating the effect of the 
local boundary layer on the vents. The growth of the boundary layer 
along the vehicle was influenced by the interaction with flanges pro- 
truding into the flow over the body and longitudinal corrugations in the 
vehicle surface. 

To determii^ these discharge coefficients, a series of tests were 
conducted in the Lewis 8x6 foot supersonic wind tunnel (ref. 1). Since 
it was impractical to test the full scale vehicle in the wind tunnel, the 
boundary layers that the vents would see in free flight were calculated, 
and these boundary layers were artifically produced in the wind tunnel 
using thickening devices. The discharge coefficients of the full scale 
vents were then determined using these artifically thickened boundary 
layers in the wind tunnel test. In addition, tests were conducted on a 
3 percent Centaur shroud model. These tests measured the pressure 
distribution on the model, especially in the vicinity of the flanges, and 
determined the effect of the flanges on the boundary layer. Calculations 
similar to those made for the full scale Centaur shroud were needed for 
the 3 percent model to estimate the boundary layer height. From these 
results the flange on the 3 percent model was scaled to the boundary 
layer thickness rather than geometrically. That is, the flanges on the 
3 percent model extended into the boundary layer the same fraction of 
the boundary layer thickness as the flanges on the full scale vehicle. 

This study reports calculations that were made to predict the bound- 
ary layer parameters along the Centaur launch vehicle and the 3 perceiit 
model in the transonic flight regime. The analysis included calculations 
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of the laminar boundary layer, estimations of the laxumar-tui'bulent transi- 
tion locations and calculations of the turbulent ;3oundary layers. In addi- 
tion, a separate analysis was required to determine the effect of the 
flange on the boundary layer, and the longibidinal corrugations in the 
vehicle surface had to be taken into account. The analysis used to cal- 
culate die boundary layer for the 3 percent model was similar to that 
used for the full scale vehicle. However, die 3 percent model did not 
have the corrugated surface. 

Calculations for both the full scale vehicle and 3 percent model 
were made for Mach numbers between 0. 8 and 1. 56. From the flight 
schedule of the vehicle the free stream static pressure at Mach 0. 8 was 
6. 1504x10^ N/m^ and at Mach 1. 56 was 2. 4318x10^ N/m^. Since the 
results of the testing and therefore the analysis were needed before the 
proof fU^t, little time was available for development of elaborate ana- 
lytical techniques. Therefore, techniques that were available or could 
be quickly devised were used in the analysis. Figures 1(a) to 1(c) show 
the dimensions of the vehicle and the dimensions and locations of the 
flange and surface corrugations. The approximate methods used here 
for the effect of the corri^tions and flanges on the boundary layer could 
be applied to similar problems in other areas. 

METHOD OF ANALYSIS 

The boundary layer properties were calculated using the finite 
difference computer program of reference 2 , nd a finite difference con: - 
puter program based on the turbulence model of reference 3. The two 
computer programs gave essentially the same results. Both the laminar 
and turlnilent boundary layer properties were calculated with these pro- 
grams usit^ the appropriate viscosity. The transition location was esti- 
mated to occur during the first adverse pressure gradient on the nose 
of the vehicle and at the shoulder of the 3 percent model, figure 1. As 
would be expected, when dealing with a body the size of the full scale 
vehicle, the results at the downstream locations near the vents were 
not particularly sensitive to the choice of any reasonable transition loca- 
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tion. On the 3 percent model the comparison of die boundary layer 
parameters with the test results served as a indirect check on the 
transition location estimation. The boundary layer analysis did not 
include a separate i^ock^-boundary layer analysis but calculated the 
properties by smoothing the surface pressure distribution throu^ the 
interaction regicm to the extent that the boundary layer program could 
be used throu^ the interaction region. 

For the areas of the Centaur shroud with the corrugated surface, 
figure 1, a method was needed to determine the effect of the surface <m 
the growth of the boundary layer. The problem is a three dimensional 
one. Becmise of the need for immediate results, however, an effort 
was made to estimate the effect of the corrugations rather than attempt 
a more rigorous three dimensional analvsis. The technique used in 
this study to estimate the effect of the corrugations was to determine 
an equivalent surface roughness that would cause the same thickening 
of the boundary layer as the corrugations. The boundary layer program 
described in reference 2 was capable of calculating the boundary layer 
properties for a surface that could be described by an equivalent sand 
rou^mess. The ratio of the actual corrugated surface area to projected 
area was determined from the geometry of the corrugations. It was 
then assumed that the skin friction would increase by approximately the 
same amount. Therefore, from the relationship between sand rou^- 
ness and local skin friction presented in reference 4 for flat plates, an 
equivalent sand roi^mess was determined and used in the boundary 
layer program. To summarize, the effect of corrugations on the bound- 
ary layer was estimated by assuming that an equivalent sand roughness 
could be determined that would approximate the effect of the corruga- 
tions. 

The effect of the flange interacting with the boundary layer had to 
be determined independently. The boundary layer program was termi- 
nated at the flange, the effect of the flange on the boundary layer was 
then calculated, and the boundary layer program restarted using the 
boundary layer properties determined from the interaction analysis. 

The method used to determine the influence of the flange was essential- 
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ly a control volume analysis. In this approach, the details of the flow 
in the interaction region are not considered. The boundary layer pro- 
perties downstream of the interaction region are determined from the 
upstream properties by usii^ conservation of mass and momentum 
across the interaction region. In this analysis, it is assumed that in 
the interaction region the loss in momentum as a result of the boundary 
layer interacting with the flange is sufficiently large that other losses 
in momentum may be neglected. Therefore, the interaction analysis 
primarily involves obtaining an estimate of the drag on the flange and, 
from this value of drag, obtain the change in momentum of the boundary 
layer through the interaction region. 

The drag coefficient for the flange was determined using a some- 
what different technique for the supersonic and subsonic cases. For 
the subsonic cases, the drag coefficient was determined from empir- 
ical results, such as those in reference 5, for rectangular bodies at 
the average Mach number seen by the flange. For cases where the 
local flow was supersonic, the drag coefficient for the flange was 
determined using the average Mach number determined not from the 
undisturbed local boundary layer but from the Mach number tliat would 
occure behind an oblique shock seeing the local edge conditions. This 
modification for the supersonic cases was used because for supersonic 
flow a separation region is produced upstream of the flange. For this 
situation the shock shape produced by the flange would be more accurate- 
ly described as resulting from a wedge than a rectangular body. The 
assumption here is that the drag behavior is similar to the well known 
reduction in drag experienced by a blunt body in supersonic flow when 
a spike is added to the nose, or in a two-dimensional case, when a thin 
protruding plate is added to a thick plate. In the analysis under con- 
sideration, the vehicle surface acts as the protruding plate. In calculat- 
ing the drag for the supersonic cases, the extent of the separation region 
upstream of the flai^e was estimated from an empirical correlation given 
in reference 6. From this information the amount of flow turning and the 
oblique shock strength were determined. The edge Mach number was 
then calculated from conditions behind the oblique shock and this Mach 
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number was used as in the subsonic cases for determining the drag 
on the flange. 

After determining the drag on the flange, the boui:dary layer down- 
stream of the interaction is determined from the incoming flow condi- 
tion and the value of the flange drag. The boundary layer calculations 
were then restarted using an essentially equilibrium boundary layer 
with the momentum defect calculated from the interaction analysis. 

No attempt was made to determine the exact details of the boundary 
layer flow in the separation and reattachment regions immediately up- 
stream and downstream of the flange. 

DISCUSSION OF RESULTS 
Vehicle Boundary Layer 

The boundary layer analysis for the full scale Centaur shroud vehi- 
cle was made at four vehicle Mach numbers, 0. 8, 0. 9, 1. 0 and 1. 56. 

The free stream static pressure and temperature vs. the vehicle Mach 
number for the vehicle flight schedule is presented in figure 2. This 
shows the decrease in static pressure and temperature with the Increase 
in speed and altitude. In figures 3(a) to (e) the pressure coefficient on 
the vehicle body is given for five Mach numbers as a function of vehicle 
length. This was obtained from a potential flow analysis with no flanges. 
From this information the input parameters are determined for calcula- 
tion of the boundary layer properties. 

Figures 4(a) to .(d) show the results of the boundary layer analysis 
and the control volume approach for analyzing the boundary layer-fiai^e 
interaction. Shown on the figures are the boundary layer displacement 
thickness vs. the longitudinal position along the vehicle. The step change 
in the displacement thickness corresponds to the location of the flange on 
the vehicle. 

Figure 5 shows the comparison of the analysis with the results of 
the wind tunnel exper*. .ents reported in reference 1. The data on this 
plot represent measurements taken on the wind tunnel wall using a full 
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scale model of the flange. The boundary layer upstream of the flange 
was generated to duplicate the full scale boundary layer predicted by 
the analysis. The downstream displacement thickness is compared at 
various freestream Mach numbers. There results show a fair agree- 
ment between the analysis and the experimental measurements. The 
largest difference in the comparison occurs at the lower Mach numbers. 
As discussed in the Method oi Analysis section, determination of the 
flange drag coefficient lor modeling the boundary layer interaction at 
the lo'wer Mach numbers may have contributed to the inaccuracies in 
the calculations. From examining the results of the Mach 1. 56 case, 
the assumption of a wedge-like interference for the flange in cases of 
supersonic flow seemed to predict the interaction correctly. 

3 Percent Model Boundary Layer 

The boundary layer analysis for the 3 percent model was made at 
four Mach numbers: 0. 8, 1.0, 1. 2 and 1. 56. The model contours and 
pressure distributions were obtained from figures 1 and 3 using the 
appropriate scaling for the model dimensions. Boundary layer calcula- 
tions were first made for the model without the rings that simulated the 
flanges on the full scale vehicle. The results of these calculations are 
shown in figures 6(a) to (d) as continuous lines. The plots show the dis- 
placement thickness as a function of longitudinal distance along the model. 
Also on the plots are the boundary layer measurements taken on the 
model in the wind tunnel tests. The agreement between the analysis and 
the data is very good at the supersonic Mach numbers but it deteriorated 
at the lower Mach numbers. This loss in accuracy at transonic Mach 
numbers could result from disturbances feeding foinvard from the solid 
rocket boosters that were located just downstream of the boattail section 
of the shroud. The boundary layer analysis would not account for this 
upstream influence. 

Also shown on figure 6 are the analytical and experimental results 
for the 3 percent model with a 0. 047 cm high ring to simulate the effect 
of the full scale vehicle flange. The 0. 047 cm height of the ring re- 
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presents the fraction of the boundary layer on the model equivalent to 
the distance the flange extends into the boundary layer on the vehicle. 

For this case the boundary layer growth up to the rii^ is the same as 
in the no ring co.ifiguration. A step change in the boundary layer thick- 
ness parameters is shown at the ring location. The results of the bound- 
ary layer calculations are then shown to the aft of the model where the 
analysis is compared with the results of the wind tunnel measurements. 
Here, as in the case without the rings in place; the comparison is better 
at the supersonic Mach numbers than at the lower Mach numbers. The 
good agreement between the analysis and measurements at the forward 
probe location indicates that the laminar-turbulent transition location 
was predicted correctly. 

An even more critical comparison between the boundary layer analy- 
sis and the experimental wind tunnel measurements can be made by com- 
paring the predicted and measured boundary layer profiles. Figures 7 
to 10 show the calculated boundary layer velocity profiles along with the 
measured values for the four Mach numbers under consideration. These 
plots show the excellent agreement at the supersonic Mach numbers. 

The difference between the analytical and measured profiles at the lower 
Mach numbers occurs in the outer portion of the boundary layer. These 
results are consistent with the suggestion that the measured profiles 
are being distorted by some transonic flow interaction on the model at 
these Mach numbers. It should be noted that there is a change in the 
scale of the ordinate between the forward boundary layer profiles and 
the aft profiles. The thickness of the aft boundary layer profiles is 
approximately twice that of the forward profiles. 

SUMMARY OF RESULTS 

To estimate the boundary layer properties over vents on the stand- 
ard Centaur shroud, boundary layer analyses were used that utilized 
approximate techniques to include the effects of flange -boundary layer 
interactions and loi^itudinal surface corrugations. Calculations were 
made for the full scale vehicle, a 3 percent wind tunnel model and a full 
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scale portion of the shroud tested in the wind tunnel. Comparisons be- 
tween wind tunnel data and analysis indicated: 

1. The results of the full scale vehicle calculations show good agree- 
ment with the boundary layer measurements downstream of the flange 
for supersonic flow but gave smaller displacement thicknesses than the 
e:q)erimental data for the subsonic free stream cases. 

2. Comparison of the 3 percent model calculations with the wind 
tunnel data show a trend similar to that of the full scale flange wind 
tunnel data, that is, good agreement at supersonic Mach numbers with 
less agreement for the aft boundary layer location for the subsonic Mach 
numbers. 

3. The overall results show that the approximate boundary layer 
techniques described here could be used in other similar problems, where 
the exact details of the boundary layer flow in the flange interaction region 
are not needed. 
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APPENDIX - SYMBOLS 
P -P«o 

Cp pressure coefficient = 

Poo^oo 

2 

M Mach number 

p static pressure 

r radius 

t static temperature 

u local flow velocity 

X axial distance 

y distance normal to surface 

p fluid density 

6* boundary layer displacement thickness 

Subscripts: 

e boundary layer edge 

r sea level 

free stream 


00 
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Figure 2 . - Free streim static pressure and temperature vs vehicle Mach number. 
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ta) Surface pressure coefficient vs distance along vehicle at - 0. 8. 




Ic) Surface pressure coefficient vs distance along vehicle at Mk • 1. 0. 

Figure 3. 
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It) Surbtt pressure coefncient vs distinct along vehicit at • 1. S6. 
Figure 3. • Concluded. 


Dl9lacmMtlhlckiMU, 6*-cni 





|i) Srintfiry 4is^lactfiiifit ttiidmtss vs tfisiHKt atong vihicft tor ^ • 1.51 

Figurtl 


ORIGINAL PAGE IS 
OF POOR QUALITY 



Figure 5i - Coinpirison of experimental and analyticat displacement 
thiduiesses doumstreiim of the flange for mtious free stream 
Mach numbers. 



(ci Boundaiy layer displacement thickness vs distance along 3 per* 



Figure 6. 



0 .4 .8 1.2 0 .4 .8 1.2 0 .4 .8 1.2 

U/Ue 

<3) Upstream profile. I>) Aft profile w/o ring. Icl Aft profile with ring. 


Figure 7. ■■ Boundary layer velocity profiles on 3 percent model for M - 0.8. 
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Figure 9. - Boundary layer velocrty profiles on 3 percent model for M«> • 1*2. 



Figure 10. - Boundary layer velocity profiles on 3 percent model for Mqo * 1. 36. 


